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Abstract 

We show that monoenergetic ion beams can be accelerated by moderate Mach number collision- 
less, electrostatic shocks propagating in a long scale-length exponentially decaying plasma profile. 
Strong plasma heating and density steepening produced by an intense laser pulse near the critical 
density can launch such shocks that propagate in the extended plasma at high velocities. The 
generation of a monoenergetic ion beam is possible due to the small and constant sheath electric 
field associated with the slowly decreasing density profile. The conditions for the acceleration of 
high-quality, energetic ion beams are identified through theory and multidimensional particle-in- 
cell simulations. The scaling of the ion energy with laser intensity shows that it is possible to 
generate ~ 200 MeV proton beams with state-of-the-art 100 TW class laser systems. 

PACS numbers: 52.38.Kd, 52.35. Tc, 52.35.Mw, 52.38.Dx, 52.65.Rr 
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The electric and magnetic fields excited in a plasma by laser pulses allow for the ac- 
celeration of ions to high energies over short distances pQ. Such accelerated beams are of 
interest for a broad range of potential applications: cancer therapy [21 E], isotope genera- 
tion [3], proton radiography [5], and fast ignition [6]. Many applications require low energy 
spread (1-10% FWHM) and low emittance ion beams. Radiotherapy applications require a 
relatively high beam energy, in the range of 100 — 300 MeV/a.m.u. [7j. 

A significant effort has been devoted in recent years to determining the optimal param- 
eters for ion acceleration. The two most studied mechanisms are target normal sheath 
acceleration (TNSA) [5] and radiation pressure acceleration (RPA) [9]. However, the pro- 
duction of high-energy and high-quality ions beams remains a challenge. A third mechanism 
for accelerating ions to high energies in a laser-produced plasma is shockwave acceleration 
(SWA) [TOTHl] . Here, a shock propagates in the plasma with a velocity v s h which can re- 
flect ions from the background plasma to a velocity v ions ~ 2v s h- Previous theoretical and 
numerical studies focused on the conditions under which shocks are formed in solid targets 
using extremely high laser intensities and the accelerated ion spectrum was broad [TTI IT2] . 
Motivated by recent experimental results on monoenergetic acceleration of protons [15J, we 
consider SWA in near critical plasma density targets at modest laser intensities. 

We first derive the conditions for ion reflection from an electrostatic collisionless shock 
as a function of the initial density and temperature profile, and for relativistic tempera- 
tures. We characterize and optimize a new acceleration regime which combines SWA with 
a special regime of TNSA. In this regime, the shock propagates through an extended ex- 
ponentially decreasing density profile, where the sheath fields are small and constant |16j . 
and the shock-accelerated ions maintain a narrow energy spread. We show, through multi- 
dimensional particle-in-cell (PIC) simulations, that shocks satisfying these conditions can be 
driven in laser-plasma interactions at near critical density and confirm our theoretical scal- 
ing of the final proton energy with laser intensity (or normalized vector potential a^). Our 
results establish the conditions to achieve the beam energy and quality required for medical 
applications in near future experiments with readily available laser systems (ao ~ 10). 

To study shock formation and ion acceleration, we consider the interaction of two 
plasma slabs (denoted by plasma 1 and plasma 0) with electron temperature ratio = 
T e i/T e0 and density ratio T = n el /n e0 (Fig. [T] a). Electrostatic shock structures can 
be generated as a result of the expansion of plasma 1 (downstream - region behind the 
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shock) into plasma (upstream - region ahead of the shock) [13j E]> with the dissipa- 
tion provided by the trapped particles behind the shock and, for strong shocks, by the 
ion reflection from the shock front [18J. The nonrelativistic theory, whereby an elec- 
trostatic shock is supported by regions/slabs of arbitrary temperature and density ra- 
tios, has been outlined in [13]. We have generalized this theoretical framework for rel- 
ativistic electron temperatures [20] and use it here in order to study the optimal con- 
ditions for ion reflection by the shock. The generalized nonlinear Sagdeev potential is 
given by ^(y?) = P;(</?,M) - P el (<p, O, T, /i e0 ) - P e0 ((p, T, /x e0 ), where P el {ip, O, V, /i e0 ) = 



re/(i + njoweV^t/We] [JT ^e-^V^ + ^eo) 2 - 1 + e-^/ Q (<rV^T - 

Log[cr + \fa 2 — 1]) — l| is the downstream electron pressure, P e0 (ip,T, fi e0 ) = 1/(1 + 
r) | (fi e0 / Ki[fi e0 \) d7e _/ie07 A/(7 + ip/ fi e o) 2 — 1 — l| is the upstream electron pressure, 
and Pj(y9,M) = M 2 (l — a/1 — 2(p/M 2 ) is the ion pressure for the assumption of cold ions 
and relativisitic Maxwellian electrons, with /z e0 = m e c 2 /kBT e0 the inverse of the normalized 
electron temperature. Here, M = v s h/c S Q is the shock Mach number, c s q = (kBT e o/rrii) l l 2 
is the upstream sound speed, ip = ecj)/kBT e0 is the electrostatic potential energy across 
the shock front normalized to the upstream thermal energy, fee is the Boltzmann constant, 
rrii and m e are the ion and electron mass, a = 1 + ip/fi e o, and K\ is the modified Bessel 
function of the second kind. In the ultra-relativistic limit, yU e0 1, P e i(ip,Q,T, fi e0 ) = 
</T> e0 (l - <p/e) + + 0]/[(l + r)6], and P e0 (<^, 6, T, /i e0 ) = <p(l - // e0 )/(l + T). Shock 
solutions can be found for ty(<p) < [19]. Ion reflection from the shock front will occur when 
the electrostatic potential across the shock exceeds the kinetic energy of the upstream ions, 
i.e. ip cr = M 2 r /2, which defines the critical Mach number 



Two-dimensional (2D) PIC simulations, performed with OSIRIS 2.0 [21], of the interac- 
tion of two semi-infinite plasma slabs with different temperature and density ratios are in 
good agreement with theory (Fig. [l]b). Ion reflection can occur for moderate Mach number 
shocks provided that r > 1 and ~ 1. At high density ratios T > 4, the expansion of 
the two slabs (initially at rest) is sufficient to form the shock and reflect the ions. At lower 
density ratios, the plasma slabs need to have an initial relative drift in order to reach M cr 
for ion reflection. 

In more realistic configurations, where finite slabs are considered, it is important to 



address the role of competing accelerating fields. As hot electrons expand into vacuum, 
TNSA fields will develop at the plasma-vacuum interface accelerating the upstream ions 
to a given velocity Vq. The shock will then reflect the upstream ions to a velocity Vi ons ~ 
2M cr c S Q + Vq. This is shown in Fig. [2] by simulating the interaction of two finite plasma 
slabs with F — 10 and = 1. For an abrupt plasma- vacuum transition, the electrostatic 
field in the sheath at the rear side of the upstream plasma introduces a chirp in vq [llj, 
broadening the ion energy spectrum as typical of TNSA [HI [16] (Fig. [2] a). This sheath 
field can be controlled by using an exponential plasma profile with scale length L g , which is 
characterized by a constant electric field at early times (t <C 4L g /c s0 ) [22] given by 

k b T e0 

&TNSA — j — , 

eL g 

as illustrated in Fig. [2] b, where we replace the low density slab of Fig. [2] a with an 
exponentially decreasing profile. It can be observed that the shock accelerated ions are able 
to cross the sheath region while preserving their narrow energy spread, thus indicating a 
configuration suitable for the generation of monoenergetic ion beams. 

The conditions for shock formation and monoenergetic ion acceleration can be obtained 
in practice from the interaction of a moderate intensity laser pulse with a tailored plasma 
density profile (see Fig. [3] a). For near critical density plasmas a significant fraction (> 20% 
[21?] 121]) of the laser energy can be absorbed maximizing electron heating, and therefore 
ion acceleration. In the relativistic regime, the electron temperature, 3/csT e = e e , can be 
estimated by equating the plasma electron energy density to the absorbed laser energy den- 
sity, 3aon c L target kBT e = r]Iri aser , where r\ is the absorption efficiency and the relativistically 
corrected critical density a n c has been used, yielding 

T e [MeV] ~ 0.026r/a - w[ps] (3) 
L target [mm\ 

For typical picosecond scale laser pulses with relativistic intensities, oq > 1, and target 
size L target < 1 mm, strong heating to MeV temperatures can occur, leading to high shock 
velocities and high reflected ion energies. 

For the production of a monoenergetic ion beam, the shock velocity should be uniform, 
which relies on uniform heating of the plasma electrons. That is achieved by allowing the 
heated electrons to recirculate in the target before shock formation [11] [25]. From [26] 
the shock formation time is ~ 47r/u P i for M ~ 1. For MeV electrons to recirculate at 



least once in the target, the target size should be limited to L target < (mj/m e ) 1 / 2 A for 
critical density targets. The shock will efficiently reflect a uniform ion population if the 
expanding ion velocity, vq = (c 2 /L 5 )t, is much smaller than the shock velocity by the time 
the shock is formed, i.e. if L g 2(c 2 /t> s / l c)(m i /m e ) 1 / 2 Ao. For a symmetric target expansion 
(Ltarget < 2L 9 ) and low Mach number shocks (M > 1), the optimal target scale length for 
uniform electron heating and ion reflection is 

A / rrii \ 1/2 

L 90 « IT — ■ (4) 



5 [MeV] ~ 2M 2 T e0 [MeV] + Mff ^ (^oW) 3/2 + 



acc 

LgQ 



2 V m, , 

Stable SWA requires a shock width (which is close to the laser spot size Wq) larger than the 
transverse expansion of the plasma, at c s , during the acceleration. Assuming an isothermal 
expansion, this condition yields W > L g0 /M cr . 

The final ion energy is given by the combination of SWA with the uniform expansion of 
the upstream plasma. The final relativistic ion velocity is Vi ons = (v' sh + v )/(l + v' sh vo/c 2 ), 
where v' sh = (2Mc s o)/(l + M 2 c 2 /c 2 ) is the velocity of the reflected ions in the upstream 
frame and vq is the upstream velocity at the shock acceleration time t acc . Taylor expanding 
Vi ons for c g o/c -C 1, the proton energy for optimal conditions is 

J mi/m e 

(5) 

To explore the proposed generation of high-quality ion beams to 10s — 100s MeV from a 
laser-driven electrostatic shock we have performed 2D OSIRIS simulations. The simulation 
box size is 3840 x 240 (c/oo ) 2 with 12288 x 768 cells, 9 — 36 particles per cell per species, 
cubic particle shapes, and current smoothing. We start by modeling the interaction of 
a Gaussian laser pulse, duration of 1885WQ -1 (FWHM) and infinite spot size with a pre- 
formed electron-proton plasma profile with a linear rise over 10Ao, and exponential fall with 
L g = 20Ao (according to Eq. Q). Increasing laser intensities (a = 2.5 — 20) have been 
used and the peak density of the plasma (n p /n c = 2.5 — 10) changed to compensate for an 
increased relativistic transparency. The laser pulse is highly absorbed as it interacts with the 
near critical density plasma (~ 60 % absorption) and stopped at the critical density surface 
causing a local steepening and leading to a density spike with 3 — 4 times the background 
density (Fig. |3]b). This density spike leads to the onset of shock formation, around t ~ 
4500 Uq 1 (530 1 after the laser interaction finished). The shock structure has a strong 



localized electric field at the shock front, with a measured thickness of L s h ~ 4Ad = 10c/a;o, 
where = ^kBT e /ATtn p e 2 is the Debye length, much smaller than the mean free path for 
particle collisions (L sh < A ei ~ c/y e , ~ 2 x 10 8 \ D , Xa ~ c a o/va ~ 2 x 10 2 \ D , for T e — 1 
MeV, Ti = 100 eV, and n e = m = 10 21 cm" 3 ). 

After the formation, the shock maintains a near uniform velocity, with a Mach number 
in the upstream reference frame M = (v s h — vq)/c s o ~ 1.7, in good agreement with the 
theoretical M cr for large T and ~ 1, M cr ~ 1.5 — 1.8 (Fig. [j]b). The temporal ion 
phase-space evolution is shown in Figures [3] c-h for the cases of a = 2.5, 10, and 20, where 
it is possible to observe the self-similarity of the interaction. The shock is able to reflect 
the cold, uniformly expanding ions from the back of the target generating a beam with an 
energy of 31, 165, and 512 MeV, respectively. The upstream ion temperature measured 
during the acceleration is relatively small (100 keV for ao = 2.5 and 1 MeV for a$ = 20). 
The uniform shock and upstream velocities lead to a reflected ion beam with a low total 
energy spread of ~ 10% (FWHM) and an average slice energy spread of 4 % (FWHM). The 
laser to ion beam energy conversion efficiency is measured to be 2 — 3 % in all simulations. 
The fraction of upstream ions reflected by the shock ranges between 10-20%. Assuming 
cylindrical symmetry, the total number of accelerated ions as inferred from the simulation 
is given by N ions ~ 10 10 (W / o[/W m ]) 2 /Ao[/ini], where Wo is the laser spot size, ideal for most 
applications. For instance, in radiotherapy ~ 10 8 ions per bunch are used in multi-shot 
treatment and ~ 10 11 ions per bunch in single shot treatment [21 E]. 

We have confirmed that our picture for SWA is still valid for a finite laser spot size by 
performing 2D simulations under the same conditions, with a$ = 2.5 and a super-Gaussian 
spot size Wo = 16Ao- A monoenergetic ion beam with 28 MeV and a narrow energy spread 
of 9 % similar to that shown in Fig. |3]d was produced. This ion beam has a small divergence 
of 4° half angle (Fig. gi). 

It is important to note that the acceleration of ions by the shock occurs after the laser 
has fully interacted with the plasma; the accelerated beam properties do not depend on the 
exact laser pulse profile and are not significantly affected by laser-plasma instabilities in the 
front of the target, such as filamentation. The beam quality depends mainly on the target 
profile. Simulations performed outside the optimal parameter range (Eq. Q) led to an ion 
bunch with a larger energy spread and/or less energy (using L g = 40 Ao resulted in a 17 MeV 
ion bunch with an energy spread of ~ 30 %). 
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The scaling of this scheme with the laser /plasma parameters was investigated by com- 
paring the simulation results with our theoretical estimates. The electron temperature is 
observed to scale linearly with the laser amplitude (Fig. [4] a), which is consistent with Eq. 
(|3j) for a laser to electron coupling efficiency rj = 0.51 (also consistent with our measured 
laser absorption). For the same target profile (with a relativistically corrected peak density), 
T and G are fixed and so is M cr , thus the proton energy will depend mainly on the electron 
temperature. This is confirmed by the measured proton energy scaling with a , which is in 
good agreement with Eq. ([5]), for an acceleration time of t acc = 5500CJQ -1 (consistent with the 
average acceleration time in our simulations (Fig. |4]b)). At low intensities the acceleration 
is dominated by shock reflection (first and second terms of Eq. ([5])), but at higher intensities 
the contribution from the ion expansion (third term of Eq. (|5])) also becomes important, 
leading to a transition from a scaling with o?J 2 to a^. This favorable scaling allows for the 
generation of high quality ~ 200 MeV proton beams, required for medical applications [7] , 
with a 100 TW class laser system (a = 10). 

In conclusion, we have presented a scheme for the generation of monoenergetic ion beams. 
Ions are accelerated by an electrostatic shock driven in an exponentially decaying plasma 
with a peak density close to critical density. The interaction of an intense laser pulse with 
such plasmas results in density steepening and strong electron heating which facilitates the 
formation and stable propagation of a moderate Mach number collisionless shockwave. The 
narrow spectrum of ions reflected from such a shock is preserved in this tailored plasma 
profile since the sheath field is constant and small. The high-quality and favorable scaling 
of the process with laser intensity pave the way for the generation of the ion beams required 
for medical applications with readily available laser systems. 
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FIG. 1: a) Schematic representation of the interaction between two plasma slabs with different 
temperature and density (red and blue), which leads to shock formation (black line represents the 
electrostatic potential) and ion reflection, b) Critical Mach number for ion reflection as a function 
of the density ratio T and temperature ratio between the two plasma slabs /regions, for T e o = 1 
keV (dashed line [13]) and T e o = 1.5 MeV (solid line Eq. Q). The symbols indicate the simulation 
values for the nonrelativistic (+) and relativistic (o) electron temperatures, obtained by measuring 
the speed of the shock structure (density jump or electrostatic field) when ion reflection is observed. 
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FIG. 2: Ion phase-space at 7700 u)~^ after the interaction of two plasma slabs with initial temper- 
ature of 1.5 MeV, = 1, and Y = 10. In a) a flat density profile is used for the low density slab, 
whereas in b) it is replaced by an exponential profile. The black lines indicate the initial plasma 
density profile and the blue lines indicate the early (t = 560 ^ p i) longitudinal electric field. The 
thin black lines indicate the integrated ion spectrum ahead of the shock. 
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FIG. 3: a) Schematic representation of SWA driven by the interaction of a laser (green) with a 
near critical density plasma (dashed). The laser heats up the plasma electrons and steepens the 
density profile at the critical density (red), driving a shock which reflects the upstream expanding 
ions, b) Transversely averaged density and field structure of the shock driven by a laser with 
arj = 2.5. The density profile is shown at t = (black), at the interaction of the peak of the laser 
with the critical density (red), and after shock formation (dashed red), t = 6560 Wq 1 , together 
with the longitudinal electric field (blue). The uniform TNSA field is in good agreement with Eq. 
([2]) for the measured T e o = 1.6 MeV (Etnsa = 0.025m e cuo/e). c - h) Ion phase-space evolution 
for ao = 2.5 (c, d), 10 (e, f), and 20 (g, h). Shock position is indicated by the arrows and the black 
lines indicate the final integrated spectrum of the reflected ions, d) Momentum distribution of the 
accelerated ion beam for a finite laser spot size Wo = 16 An. 
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FIG. 4: Scaling of a) electron temperature and b) shock accelerated proton energy with laser 
ao- The electron energy distribution obtained in OSIRIS is well fitted by a relativistic Maxwellian. 
The obtained scalings are consistent with Eq. ^ for r/ = 0.51 and Eq. ^ for t acc = 
respectively. 
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